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The problem of anomalous diffusion in momentum (velocity) space is considered 
based on the master equation and the appropriate probability transition function 
(PTF). The approach recently developed for coordinate space in [1], is applied with 
necessary modifications to velocity space. A new general equation for the time evo- 
lution of the momentum distribution function in momentum space is derived. This 
allows the solution of various problems of anomalous transport when the probability 
transition function (PTF) has a long tail in momentum space. 

For the opposite cases of the PTF rapidly decreasing as a function of transfer 
momenta (when large transfer momenta are strongly suppressed), the developed ap- 
proach allows us to consider strongly non-equilibrium cases of the system evolution. 
The stationary and non-stationary solutions are studied. 

As an example, the particular case of the Boltzmann-type PT- function for 
collisions of heavy and light particles with the determined (prescribed) distribution 
function, which can be strongly non-equilibrium, is considered within the proposed 
general approach. The appropriate diffusion and friction coefficients are found. 
The Einstein relation between the friction and diffusion coefficients is shown to be 
violated in these cases. 
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I. INTRODUCTION 

Interest in the anomalous diffusion is prompted by a large variety of applications, i.e., 
semiconductors, polymers, some granular systems, plasmas under specific conditions, various 
objects in biological systems, physical-chemical systems, and others. 

Deviation from the linear in time dependence < r^(t) >~ t of the mean-squared displace- 
ment has been experimentally observed, in particular, under essentially non-equilibrium con- 
ditions and in some disordered systems. The mean-squared separation of a pair of particles 
passively moving in a turbulent flow anomalously increases (according to the Richardson's 
law) with the third power of time [2]. For diffusion typical of glasses and related complex 
systems [3], the observed anomalous time dependence is slower than linear. These two types 
of anomalous diffusion are referred to as superdiffusion, < r^{t) >~ (a > 1), and subdif- 
fusion (a < 1), [4]. To describe these two diffusion regimes, a number of efficient models and 
methods have been proposed. The continuous time random walk (CTRW) model of Scher 
and MontroU [5] , leading to strongly subdiffusive behavior, provides a basis for understand- 
ing photoconductivity in strongly disordered and glassy semiconductors. The Levy-flight 
model [6], leading to superdiffusion, describes various phenomena such as self-diffusion in 
micelle systems [7] , reaction and transport in polymer systems [8] , and has been even to the 
stochastic description of financial market indices [9]. For both cases, the so-called fractional 
differential equations in coordinate and time spaces were successfully applied [10]. 

However, recently a more general approach has been proposed in [1,11], which reproduces 
the results of the standard fractional differentiation method (when it is applicable) and allows 
the describe of more complicated cases of anomalous diffusion. In [12], this approach has 
also been applied to diffusion in a time-dependent external field. 

In this paper, the problem of anomalous diffusion in momentum (velocity) space is con- 
sidered. In spite of the formal similarity, diffusion in momentum space is very different 
physically from coordinate space diffusion. This is obvious because momentum conservation 
which takes place in momentum space has no analogy in coordinate space. 

Some aspects of anomalous diffusion in velocity space have been considered in several 
papers [13-18]. In general, in comparison with anomalous diffusion in coordinate space, 
anomalous diffusion in velocity space is poorly studied. To our knowledge, there still is no 
corresponding way to describe anomalous diffusion in velocity space self-consistently . 
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In this paper, a new kinetic equation for anomalous diffusion in velocity space is derived 
(see also [19,20]) based on the appropriate expansion of the PTF (in spirit of the approach 
proposed in [1] for the diffusion in coordinate space) and some particular problems are 
studied on this basis. 

The paper is organized as follows. The diffusion in velocity space for the cases of the 
normal and anomalous behavior of the PT-function is presented in Section II. Starting from 
argumentation based on the Boltzmann-type of the PTF, we derive a new kinetic equation 
which then can be applied not only to the Boltzmann-type processes, but also to a wide 
class of processes occurring in nature with other types of PT-functions. The particular cases 
of anomalous diffusion for hard spheres collisions with the specific power-type prescribed 
distribution function of light particles are analyzed in Section III. More general examples 
of anomalous diffusion are also considered. The non-stationary solution for the distribution 
function identical at t = to the initial distribution is found. In Section IV, the Boltzmann- 
type equation is used to consider the effect light particle drift on the PT-function and on the 
kinetic equation structure. In Section V, a short review of anomalous diffusion in coordinate 
space is presented. The method and results of this paper can be easily extended to describe 
plasma-like systems and anomalous processes in energy space. 



II. DIFFUSION IN VELOCITY SPACE ON THE BASIS OF THE 
MASTER-TYPE EQUATION 

Let us now consider the main problem formulated in the introduction, namely, the de- 
scription of diffusion in velocity space (V-space). Such a description would be based on 
a corresponding master equation describing the balance of grains coming at and from the 
point p at the instant t. The structure of this equation 

^^^^ = y"^qmq,P + q)/,(p + q,t)-H^(q,p)/.(p,t)}. (1) 

is formally similar to the master equation Eq. ( 1621] (see below) in coordinate space. Surely, 
for coordinate space, there is no conservation law similar to that in momentum space. The 
probability transition function W{q^, p') (the transition rate) describes the probability that 
a grain with momentum p' (point p') passes from this point p' to the point p per unit time 
by transferring the momentum q = p' — p to the surrounding medium. Assuming initially 
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that the characteristic transfer momentum q is smaller than p, one may expand Eq. ([T]) in 
on q to the second order. Then, we arrive at the usual form of the Fokker-Planck equation 
for the density distribution function /^(p, t) (see, e.g., [21]) 



The coefficients A^, and Ba/3 describe the friction force and diffusion, respectively. The 
indices a and (3 correspond to the coordinate axes x, y, z (for dimension s = 3). 

Since the velocity of heavy particles is low, the p-dependence of the PTF can be neglected 
for calculating the diffusion term, which in this case is a constant Bap = SapB. The constant 
B is equal to the integral 



where s is the dimension of velocity space (V-spaceO. Neglecting the p-dependence of the 
PT-function, we arrive at the coefficient Aa = (while the diffusion coefficient is constant). 

As is well known this neglect is wrong and the coefficient A^ for the Fokker-Planck equa- 
tion can be determined using the assumption that argument that the stationary distribution 
function is Maxwellian. Then, we arrive at the standard relation between the coefficients 
MTAa{p) = PaB (here M is the particle mass and T is temperature of the particles in the 
equilibrium state). This relation is an analogue of the Einstein relation in coordinate space. 
However, this argumentation is not applicable to systems far from equilibrium,. 

To find the coefficients in the kinetic equation, which are appropriate for more general 
non-equilibrium situations, e.g., for slowly decreasing PT-functions or (and) for strongly non- 
equilibrium kernels W, let us use a more general way based on a certain small parameter 
(in the simplest case, e.g., on the difference of velocities of light and heavy particles for 
Boltzmann-type collisions). To calculate the function A^ we have take into account that 
the function Vr(q, p) is a scalar and can depend only on variables (?,q-p,P- Expanding 
W{c{, p) in the product q ■ p, one arrives at the approximate representation of the functions 
W{ci, p) and W{ci, p + q). 






(3) 




(4) 





W^(q, p + q) ^ W{q) + W'{q) (q ■ p) + ^W"iq){q ■ + q'W'{q), 



(6) 
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Here we introduced the derivatives W'{q) = dW{q,q - p,p)/9(qp) |q.p=o,p=o and W"{q) = 
d'^W{q,q-p,p)/d{qp)^ |qp=o,p=o- 

Then, with necessary accuracy (which corresponds to the derivation of the well known 
Fokker-Planck equation) for the coefficient Aa we find 

^a(p) = J d'qqaqm^\q) = Pa J d'qqaqaW'{q) = ^J d'qq^W\q) (7) 

If the equality W'{q) = W{q)/2MT is fulfilled for the function W^(q, p), we arrive at the 
usual Einstein relation for the coefficients Aa and Ba/s 

MTAa{p)=Po.B (8) 

Let us check this relation for Boltzmann collisions which are described by the PT-function 
W{q,p) = WB{q,p) [11]: 

,2 



2n da 

^^q Jq/2n do 



arccos(l - ),u 



^{u' + v'-q■^r/fi), (9) 

where (p = Mv), the quantities da{X:'^)/do , /i and fb are, respectively, the scattering dif- 
ferential cross- sect ion, the mass and distribution function for light particles. The arguments 
X and u in the cross-section da /do are the angle of the particle scattering in the center- 
of- mass coordinate system and the velocity of light particles before collision. In Eq. (jO]), 
we took into account the approximate equalities for light and heavy particle scattering, 
q^ = (Ap)^ = p'^(l — cosO) and 6 — x (where p' = fiu is the light particle momentum before 
collision). For the equilibrium Maxwellian distribution /°, the equality W'{q) = W{q)/2MT 
is evident and we arrive at the ordinary Fokker-Planck equation in velocity space with the 
constant diffusion D = B/M"^ and friction /3 = B/MT = DM/T coefficients which satisfy 
the Einstein relation. 

The problem of the determination of the W^-function and the respective Fokker-Planck 
equation for the situations close to equilibrium are discussed in detail in the review [22]. 
However, even for quasi-equilibrium regimes, where long tails of the PTF-functions are 
absent, the consideration in [22] is restricted to the hard sphere interaction. The non- 
equilibrium forms of the VT-function and the drastic changes in the Fokker-Planck equation 
structure for these situations were not considered, as well as the specific situations, when 
the long tail of the I^-function exists. Therefore, anomalous diffusion in [22] is absent. As 
is easy to see, in the case of the hard-sphere cross-section and the equilibrium Maxwellian 
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distribution, fb Eq. ([9]) leads to the same result for the Wb function as has been discussed 
in [22]. 

For some non-equilibrium (stationary or non- stationary) states, the PTF can have a long 
tail as a function of q. In this case, we have to derive a generalization of the Fokker-Planck 
equation in spirit of the consideration [1,11] for coordinate case. The necessity of this 
derivation arises since the diffusion and friction coefficients in the form Eqs. dlD,© diverge 
in the limit of large q, if the kernels of the functions W, W exhibit the asymptotic behavior, 
W{q) ~ 1/g" with a < s + 2 and (or) W'{q) ~ 1/g^ with /3 <s + 2. 

Let us substitute the expansions for W (as an example, we choose s = 3, the arbitrary s 
can be considered by the similar way) into Eq. ([T]). With an accuracy up to (q ■ p)^, we find 

dfg{p,t) _ 



dt 



c^q{/.(P + q, t) [W{q) + W\q) (q ■ p) + 



\w'\q){q ■ vf + q'W\q)] - /,(p, t)[W{q) + W\q) (q • p) + \w"{q)[q ■ p)^]} (10) 
After the Fourier-transformation fg{s,t) = J j^^exp{ipr)fg{p,t), Eq. ffTOj) reads 



dt 



dSr 



(11) 



dSadSi3 

It should be noted that the variable Sq, arisen in Eq. ( ITTl) has the dimension p^;^ and is a 
formal variable for the Fourier-transformation. 
The coefficients in Eq. f lTTj) are given by 

sin (qs) 



A{s) 



dq[exp{—i{qs)) — l]iy(g) = An j dqq^ 





qs 



Wiq) 



(12) 



B^{s) = s^Bis); B{s) = -- / dqqs[exp(-z(qs)) - l]W'iq) 



Jo 



cos (qs) — 



sin{qs) 
qs 



W'iq) 



(13) 



Cafi{s) = Sa,Si3C{s) = I dqqaq/3[exp{-i{cis)) - l]W"{q) 



(14) 



C{s) = / dqiqs)^[exp{-t{qs)) - l]W"{q) 



^ / dqq'' 



2sin{qs) 2cos (qs) sin{qs) 1 

3 



q^s^ 



qs 



W'\q) 



(15) 
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For the isotropic function /(s) = f{s), one can rewrite Eq. (fTT]) as 



dfg{s,t) 



A{s)fg{s,t) + B{s)s 



+ C{s)s' 



(16) 



In the case where the PT-function W{q) and the functions W'{q) and W"{q) strongly 
decrease for large values of q, the exponents under the integrals in the functions A{s), B{s), 
and C{s) can be expanded as 



Ais)--j I dqq'Wiq); 5(s) 



dciq'^W'iq); C{s) - 0. 



(17) 



Then, the simplified kinetic equation in velocity space PT-function (which is non- 
equilibrium in the general case) is written as 



dfg{s,t) 



dfgis) 



dt "''^ ' ds ' 

where Aq = -1/6 / dqq'^W{q) and 5o = -1/3/ dqq^W'{q) is defined by Eq. 



(18) 



For C{s) = the stationary solution of Eq. ( IT6l) is given by 



fg{s) = Cexp 



ds 



s'Bis'] 



Cexp 



A^ 
'2Bo 



(19) 



The respective normalized stationary momentum distribution is written as 



NgB'f 



■exp I 



(20) 



Therefore, the constant in Eq. f|T9|l is C = Ng, where Ng is the density of heavy particles. 
Equation ( ITSl) and the obtained distribution are the generalization (for non-equilibrium sit- 
uations) of the standard Fokker-Planck equation to the case of normal diffusion in velocity 
space (see,e.g., [21]). The characteristic feature of these physical situations is the existence 
of the fixed (prescribed) kernel W{q,p) which is defined, e.g., by a certain non-Maxwellian 
distribution of small particles fb- To show this in other way, let us make the Fourier trans- 
formation of Eq. (fTTj) with C = and respective coefficients A and B^, 

dfg{p,t) ^ ^^^ gVg(p,t) _ ^^ d{pafg{p,t)) 
dt dp'^ dpa 

Therefore, we arrive at the Fokker-Planck-type equation with the friction coefficient /3 = —-Bo 



and diffusion coefficient D 
the Einstein relation. 



-Aq/M"^. In general, these coefficients (Eq. ( ITTl) ) do not satisfy 
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In the case of equilibrium (e.g., fb = /°, see above) for the IV-function the equahty 
W'{q) = W{q) /2MTb is fulfilled. Then, with the necessary accuracy of the order /i/M, we 
find A{s) / sB{s) = Aq/Bq = MTb. In this case, the Einstein relation between the diffusion 
and friction coefficients D = (5T /M exists, and the standard Fokker-Planck equation is valid. 

In fact, the approximation C{s) ^0 is practically always applicable due to the small pa- 
rameter (e.g., /i/M for the Boltzmann-type PTF, see below. Sec. III). Therefore, the general 
kinetic equation f|T6l) for the Fourier-transform of the velocity distribution function takes 
the form 

*M . Ai.)S,i.) + B^i.fJ^, (22) 
where the coefficients A{s) and -Bq(s) are determined by f|T2l) . f|T3|) . respectively. 



III. THE MODELS OF ANOMALOUS DIFFUSION IN V - SPACE 

Now we can calculate the coefficients for the models of anomalous diffusion. We will also 
estimate for the function W" and the tensor Ca^p] although, at the end we neglect these 
terms due to the small transfer momentum in the collision process. 

At first, we calculate the simple model, i.e., the system of hard spheres with different 
masses m and M ^ m, da /do = Let us suppose that light particles in the model under 

consideration are described by the prescribed stationary distribution fb = Ubipb/u^ (where 
Ubi (pbi and Mo are, respectively, the density, non-dimensional distribution, and characteristic 
velocity of light particles). The integration variable is defined by the equality ^ = {v^ + v'^ — 

q- v//i)/Mo, 

WM,v) = ^4^r d^-MO- (23) 

If the distribution (pbiO = V^'^ (7 > 1) has a long tail, we get 

W(av)- J!2_|oo _ ed'" .241 

(1 - 7) (7 - 1) 

where ^0 = {(f /^IJ^'^ + — Q ■ v//i)/Mo. 

For the case p = 0, ^0 ~^ ^0 = and we arrive at the expression for the anomalous 

PT-function W = Wa, 

WJq, P = 0) = o-^t:^ = — — -. (25) 

«V4'P J 23-27(^-l);,4-27y3-27^27-l g27-l ^ ^ 
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To determine the transport process the structure and the kinetic equation in velocity 
space, one should also find the functions W'{q) and W"{q). 

To find the functions W'{q, p) and W"{q, p) for p 7^ 0, we use full-value ^0 = 
{q'^/Afi^ + p^/M"^ — q ■ p/Mfi)/uQ and the derivatives of these functions on q • p at p = 0, 
^'0 = -1/Mfiul, and ^'^ = 0. Then 



iy'(q,p) 



W"iq,p) 



Therefore, for p = (^0 — ^ ^0) we obtain the functions 



-7-1 



(26) 



W\q) 



W"{q) 



2M/i3M|]g27+l ' 2M2/i%5^27+3 

The function A{s), according to Eq. f|T2|) . is given by 



A{s) = Atx 



POO 

/ dqq^ 


sin (qs) ^ 


Jo 


qs 



W{q) = AixCa 



dq- 



,27-3 



sin{qs) 
qs 



(27) 



(2^ 



Comparing the reduced equation (see below) in velocity space with diffusion in coordinate 
space (see the Appendix), we can establish the correspondence (27 — 1 ^ a and W{q) = 
C/q'^^^^). This means that the convergence of the integral on the right side of Eq. (!28ll 
(3d case) is provided if 3 < 27 — 1 < 5 or 2 < 7 < 3. The inequality 7 < 3 provides the 
convergence for small values of the variable q (g — )■ 0), and the inequality 7 > 2 provides 
the convergence for g — )■ 00. 

We now establish the conditions of the convergence of the integrals for B{s) and C{s) as 



Bis) = - 
s^ 



cos (qs) 



sin{qs) 
qs 



W'{q) 



(29) 



The convergence of the function B{s) exists for small values of g if 7 < 2 and for large 
values of g (g — )■ 00 for 7 > 1/2. 

Finally, the convergence for the function C{s) is defined by the equalities 7 < 2 for small 
values of g and 7 > 1 for large values of g 



C{s) 



271 



dqq^ 



2sin{qs) 2cos (qs) sin{qs) ^ 1 



W"{q) 



(30) 



q-'S'^ q^s^ qs 

Therefore, in the case under consideration, the convergence of the functions A, B, and C 
for large large values of g is defined only by the convergence of the function A, which means 
7 > 2. To provide convergence for small values of g, it is sufficient to provide convergence 
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for the functions B{s) and C, which means 7 < 2. Therefore, for the purely power behavior 
of the function fb{0^ simultaneous convergence of the coefficients A, B and C cannot 
be provided. However, in real physical models, the convergence of the coefficients in the 
integration region g — is always provided, e.g., by a finite value of v (see Eq. f l24l) ) or 
due to the non-power behavior of the PT-function W for small values of q (compare with 
the examples of anomalous diffusion in coordinate space [1]). Therefore, in the model under 
consideration, the "anomalous diffusion in velocity space" for the power behavior of W{q), 
W'{q) and W"{q) at large values of the variable q exists for the asymptotic behavior of the 
PT-function W{q — > 00) ~ l/q^'^~^ if the inequality 7 > 2 is fulfilled. In this case the 
integrals for the coefficients A{s),B{s) and C{s) converge. At the same time, the expansion 
of the exponential function in Eqs. ( fT2l) -( |T5l) under the integrals, which leads to the Fokker- 
Planck type kinetic equation, is invalid for the power-type kernels iy(q, p). 

Now let us consider the more general model for which we will not connect the functions 
W{q), W'{q), and W"{q) with a concrete form of Vr(q, p) which is in general unknown. 
In this case one can suggest that the functions have the independent one from another 
power-type g-dependence. 

As an example, this dependence can be taken as the power-type one for three functions 
W{q) = a/q°', W'{q) = h/q^, and W"{q) = cjcp , where a, /3 and r] are independent and 
positive. Then, as follows from the above consideration, the convergence of the function W 
exists if 5 > a > 3 (for asymptotically small and large g, respectively). For the function 
W\q)^ the convergence condition for asymptotically small and large values of q is provided 
if 5 > /3 > 2, respectively. Finally, for the function ly'(g), the convergence condition is 
7 > > 5 (for asymptotically small and large values of g, respectively). 

For this example, the kinetic equation Eq. ffTTj) reads 



dfg{s,t) 



/(s,t), 



(31) 



dt 



^ dsidsj 



where 




(32) 




(33) 
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fOO 

P2 = Attc / dCC^"" 
Jo 



sinC, cosC, sin( ^ 1 



(34) 



2C 6_ 

As is easy to show, Pq = ~|-Po| sgna < and Pi = —\Pi \ sgnb. Taking into account the 
isotropy in s-space we can rewrite Eq. fl3T|) in the form 

= PoS--'f{s, t) + s^-'Pi^/is, t) + s'-'P2^f{s, t), (35) 

Naturally, Eqs. f l3ip .f l35l) can be formally rewritten in momentum (or in velocity) space in 
terms of fractional derivatives of various orders. Therefore, as is easy to see, for the purely 
power behavior of the functions W{q), W'{q) and W"{q) the solution with the convergent 
coefficients exists for the values of the powers in the intervals mentioned above. In the case 
under consideration the universal type of the anomalous diffusion in velocity space exists 
if5>a>3, 5>/3>2 and 7 > rj > 5. This takes place even in the cases, when the 
functions W{q), W'{q) and W"{q) have not cut-off for small values of q. Surely, the general 
description is also valid for more complicated functions W, W, and W", characterized by 
the non-power g-dependence for small values of q. 

Now let us take into account an important circumstance; in general, there is a small 
parameter /x/M in the problem under consideration , which can simplify the description of 
velocity diffusion. As is easy to see, e.g., based on the particular cases (see Eq. ( |271) ) for 
the convergent kernels of anomalous transport, the term with the second derivative W" in 
general equations f|TTl) . f|T6|) for the distribution fg{p) is small (in comparison with the term 
with the second derivative. The similar smallness is observed for the case of normal diffusion 
in velocity space. This smallness is of the order of the small ratio fi/M of particle masses. 

Therefore, for the most physically important kernels describing anomalous velocity dif- 
fusion, the term with the second space derivative can be omitted and the non- stationary 
general diffusion equation is given by 

*M = AW/,(s.<) + B„W^?^ (36) 

or, for the isotropic case, 

*M = ^W/,(M) + BW.^ (37) 

In the case of the purely power behavior, W{q) = a/q°' and W'{q) = b/q^, we have, as above, 
A{s) = Pqs""^ and sB{s) = Pis^"^ (with the inequalities 5 > a > 3 and 5 > /3 > 2). The 
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stationary solution of Eq. (137|) (see, also ( !T9|l ) for the case under consideration is written as 



Cexp 



- / ds- 



, A{s') 
s'Bis'] 



Cexp 



Pi(a-/3 + 2) 



(38) 



where 5>a;-/3 + 2>0. 

To find the solution in the isotropic non- stationary case, Eq. (l37Il should be written as 



(39) 



where X(s,t) = lnfg{s,t). The general non-stationary solution of this equation can be 
written as the sum of the general solution of the homogeneous equation Y{s,t) (Eq. (!39|l . 
where the function A{s) is taken zero), 

1 



r(s,t) = $(0, 



ds'- 



(40) 



Here $ is the arbitrary function of the variable ^. The particular solution Z[s,t) of the 
inhomogeneous equation Eq. ( 15^ reads 



^(^,i) = /f(5 



s'B{s') 



Therefore, we find 



fg{s, t) = exp [X{s, t)] ^ exp [Y + Z] = L(0/f (s 



(41) 



(42) 



where L{C,) is the arbitrary function of ^, which should be found from the initial condition 

fg(s,t = 0) = Ms). 

The variable ^ = $,{s,t) equals 



^(s,t) = t+ ds' 



1 



,5-/3 



t + 



(43) 



so 



s'B(s') " ' Pi(5-/3) 
where c is the arbitrary constant which can be omitted due to the presence of the arbitrary 
function L. 

The general non-stationary solution for the case under consideration reads 



fg{S,t)=L(t + 



,5-/3 



Pi(5-/3) 



exp 



a-/3+2 



Pi(a-/3 + 2) 



(44) 



The inequalities for the combinations of the above coefficients are 5 > a — 13 + 2 > 0, 
3 > 5 — /3 > 0. The unknown function L can be found from Eq. (jH]) and the initial 



condition fg{s, 0) 



[s : 
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The function 4>g{s) = J pexp{ips) fg{p,t = 0) is the Fourier-component of the initial 
distribution in momentum space. Using the notation ( = s^~'^/[Pi(5 — /?)] (which means 

s(C) = [Pi (5 - (3)CY/^^-^^), we find 



L{0 = <PMO]exp 
Therefore, the time-dependent solution is 



fgis,t) = <PMC + t)] exp { ^t'^^^^T~T I exp 



Pi(a-/3 + 2) 



(46) 



(47) 



Pi(a-/3 + 2) 

where we have to express the value s{C + t) = [Pi(5 -/3)(C + t)]^/^^~^^ as the evident function 
of variables s,t: 



s(C + t) = [Pi(5-/3)(C + t)] 



1/(5-/3) 



/(5-/3) 



(48) 



or, finally. 



fg{s, t) = <Pg {[s"-^ + Pi(5 - /3) t]i/(^-«) exp 



[S^-P + Pi (5 - /3) t](«-/5+2)/{5-/3) _ 



Pi {a -13 + 2) 

It should be emphasized that the real solution for the fractional powers 1/(5-/3) and (or) 
(a — /3 + 2)/(5 — /3) exists only if Pi > 0. The limit t — ?■ oo for the solution can be identical 
to the stationary solution only for specific initial conditions. For these cases, the stationary 
solution can be non-equilibrium. 

For the power dependence of the functions W{q) and W'{q), the equation (136|) can be 
formally written in terms of fractional derivatives. 



dt 



PoD^fg{p,t)-P,{3 + ^)Dyg{p,t) + P,p^Dl^'fg{p,t) 



(50) 



where z/ = a — 3, 7 = /? — 5 and D'^'^^fg{p, t) = i f d^ s exp{—ips))saS'^ fg{s, t). 

The particular case of the collision term for anomalous diffusion in velocity space, con- 
sidered in [17], is derived in Appendix 11 on the basis of the general equations (1361) . (1471) . 



IV. DIFFUSION MODEL BASED ON BOLTZMANN COLLISIONS WITH 

DRIFT 



Let us consider the simplest case of the non-equilibrium but stationary distribution /b, 
i.e., the shifted velocity distribution. 
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The evident generalization of the PT-function w^(q, p) for this case (characterizing by 
the drift velocity Ud in the distribution function /{,) is given by 



where p = Mv. Again, as in Section 11, to find the coefficients in the kinetic equation, 
let us use the difference between the velocities of light and heavy particles. At the same 
time, the drift velocity Ud, generally speaking, is not small in comparison with the current 
characteristic velocities u and of small particles. 

To calculate the function Aa, we have take into account that the scalar function 
w^(q, p,Ud) can in general be written in the form iy(q, p,Ud) = W{q,p,Ud,l,C,,ri) (here 
/ = (q ■ pd)) and expanded on ^ = (q ■ p) and r] = (Mud ■ p) = (pd ■ p)- In fact, it is the 
expansion in the velocity v = p/M, which is small in comparison with other characteristic 
velocities q/fi, u, and, in the general case, Ud- As shown above (for the case = 0), to ar- 
rive at the simple and solvable equation for the distribution function fg, taking into account 
the smallness of v in comparison with characteristic velocities and f ^ in comparison with 
(q-v/yu)), we have approximate the function W{q,p,Uci,l,^,ri) ~ W{q,Uci,l,^,ri), because 
we are interested mainly in high values of q for anomalous transport. At the same time, after 
this approximation and expansion on and r] (for Ud 0), we can arrive (for the special case 
of kernels with the purely power-type g-dependence for small q) at the divergence of some 
coefficients of the diffusion equation. A similar situation takes place also for anomalous diffu- 
sion in coordinate space. This divergence is absent for the realistic PT-functions, which pro- 
vide the cut-off for small values of q. This cut-off has the physical nature and is not related to 
the approximation W"(q,p,Ud) = W{q,p,Ud,l,^,'n) - W{q,v = 0,Ud,l,^,r])W{q,Ud,l,^,r]). 

Let us expand the PT-function W{q,p, Ud, I, ^, r]) in ^ and t] similarly to the more simple 
case ([5]), 



W{q, p, Ud) = W{q, p, Ud, I, ^, v) - Wo{q, p, Ud, I) + dW/d^ 1^,^=0 ^ + dW/dr] 1^,^=0 V + 



Here Wo{q,p,Ud,l) = W{q,p,Ud,l,^,ri) |^,r)=o- Then, introducing the functions 



Vi{q,p,UdJ) = dW/d^ \^^r,=o,Ui{q,p,UdJ) = dW/dr] \^^^=o,V2{q,p,UdJ) = Id^W/d^'^ |g,^=o 



w^B(q,P,Ud) 




(51) 



^d'W/de k,=o e + \d^W/dr,^ |^,,=o + d^W/d^dr^ |g,,=o ^^,(52) 
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,U2{q,p,Ud,l) = \d^W/drf [5,^=0, ^^2(9,^, m^, /) = d'^W/d^dr] \^^r,=o we can rewrite Eq. §^ 
in the form 

W{q,p,Ud,l,C,r]) ~ Wo{q, p,Ud, I) + Vi{q-p) + 
UiiPd ■ p) + ^(q ■ p)' + U2{pd ■ pf + W2{q ■ p)(pd ■ p) ^ 
Wo{q,p = 0,Ud,l) + Vi{q,p = 0,Md,/)(q- p) + Ui{q,p = 0,Ud,l){pd ■ p), (53) 



Finally, we set p = in the coefficients of Eq. (1531) and omit the terms with the second order 
derivatives due to the existence of the small parameter (e.g., fi/M). This means only the 
terms of the order of ^/JlJM are essential in the expansion of W^(q, p, Ud). Let us calculate 
^(Q) P + taking into account the difference between the values of the characteristic 

momenta q, p + q, Ud or the characteristic values of velocities. 

The appropriate expansion for the function (q, p + q, Ud) with the necessary accuracy 
is written as 

1 

iy(q, p + q, Ud) ~ Vr(q, p, Ud) + {qad/dpa + i^q^q^ }W{q, p, Ud) ~ 

Z OPaOpp 

Woiq^p = 0,Ud,l) + Vi{q,p = 0,Md,/)(q ■ p) + Ui{q,p = 0,Md,/)(pd ■ p) + 

Vi{q,p = 0, Ud, l)q^ + Ui{q,p = 0, Ud, I) (q ■ Pd) (54) 

or 

W{q, p + q, Ud) ~ Woiq, p = 0, Ud, I) + 
Vi {q, p = 0, Ud, I) [q ■ (p + q)] + Ui{q,p = 0, Ud, I) [pd ■ (p + q)] • (55) 

Then the collision term of the kinetic equation can be given by the formula 

^^^^ = J dq{[Wo{q,p = 0, Ud, I) + V,{q,p = 0, n,, /)(q • p) + U,{q,p = 0, Ud, /)(pd ■ p) 

+Vi{q,p = 0,Ud,l)q^ + Ui{q,p = 0,Ud,l){q ■ Pd)] 
fg{p + q, t) - [Wo{q, p = 0, Ud, I) + Vi{q, p = 0, Ud, I) (q ■ p) + 

Ui{q,p = 0,Ud,l){pd ■ p)]/3(p,t)}56) 

After the Fourier-transformation ffTOl) . the distribution function fg{s) = 
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/ (^3ea;p(ips)/g(p, t) can be written in the form 

df (s t) f d 

= / dq{exp{-iqs)[Wo{q,p = 0,Ud,l) -iVi{q,p = 0,Ud,l){q- - 

d 

iUi{q,p = 0,Ud,l){pd ■ ■^)]fg{s,t) 

d d 
-[Wo{q,p = 0,Ud,l) -iVi{q,p = 0,Md,/)(q- — ) - iUi{q,p = 0,Ud,l){pd ■ ■^)]fg{s,t)} 

(57) 

Therefore, we arrive at the equation similar to Eq. f l5B]) . but with the evidently non-isotropic 
structure, 

= Ad{s, pd)/,(s, t) + B,(s, Pd)^/,(s, t), (58) 
where the coefficients are expressed as the integrals 

^d(s,Pd) = j dq[exp{-iqs) - l]Wo{q,p = 0,Ud,l), (59) 

Bd(s,pd) = J dq[exp{~iqs) - l]{Vi{q,p = 0,Ud,l)q + Ui{q,p = 0,Ud,l)Pd} = 

sE^(s,pd) + Pd5^'(s,Pd). (60) 

This equation will be analyzed in detail in the separate paper. 



V. CONCLUSIONS 

In this paper, the problem of anomalous diffusion in momentum (velocity) space has been 
consistently considered. The new kinetic equation for anomalous diffusion in velocity space 
has been derived for the general case of the non-equilibrium probability transition function. 

It should be emphasized that these types of anomalous diffusion are not related to the 
time dispersion (the non-local in time PT-function). Anomalous diffusion in velocity space, 
created by the non-local in time master equation has been considered in a few papers (see, 
e.g., [15], [16], [23], [24]) on the basis of the well-known Fokker-Planck kernel in p-space. 

The similar problem with the non-local in time PT-functions, which at the same time 
have long g-tails, has been considered for diffusion in velocity space separately, based on the 
results obtained in this paper. 
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The evolution of the distribution function can lead to the stationary non-equilibrium 
limit or be asymptotically time-dependent. The model of anomalous diffusion in velocity 
space is described based on the respective expansion of the kernel in the master equation. 

The conditions of the convergence of the coefficients of the new kinetic equation have been 
found for particular cases. The wide variety of anomalous processes in velocity space exists, 
since two (C is usually small) different coefficients enter the general diffusion equation even 
in the isotropic case. 

The examples of anomalous diffusion for heavy particles, based on the Boltzmann ker- 
nel with the prescribed equilibrium and non-equilibrium distribution function for the light 
particles have been studied. The PT-function for the Boltzmann kernel depends on the ap- 
propriate cross-section for scattering of two types of colliding particles and on the prescribed 
velocity distribution for one sort (light in the case under consideration) particles. The par- 
ticular case of the hard sphere interaction has been considered. The Einstein relation is 
applicable only to the equilibrium PT-function (hence, the equilibrium prescribed distribu- 
tion of hght particles). In the general case, even for the stationary solution, the Einstein 
relation between the diffusion and friction coefficients is absent. For normal diffusion, the 
friction and diffusion coefficients have been explicitly found for the non-equilibrium case. 
The stationary and non-stationary (with an initial condition) solutions of the general ki- 
netic equation have been found. In the equilibrium case, the known Fokker-Planck equation 
is reproduced as a particular case. We also have shortly reviewed anomalous diffusion in 
coordinate space. 

The kinetic equation obtained in this paper is applicable to various systems in physics 
(including plasma physics [25]), chemistry, and biology. 

VI. APPENDIX 

Appendix I. Diffusion in the coordinate space on the basis of the master-type 
equation 

Let us consider diffusion in coordinate space on the basis of the master equation which 
describes the balance of grains coming in and out the point r at the moment t. The structure 
of this equation (see, e.g., [1,11]) is formally similar to the master equation in momentum 
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space, which has been derived above, 

^^M = I dr'{iy(r,r')/,(r',t)-W^(r',r)/,(r,t)}. (61) 

Surely, for coordinate space, there is no conservation law similar to that in momentum space. 
The probability transition W{r, r') describes the probability that a grain transfers from the 
point r' to the point r per unit time. We can rewrite this equation in the coordinates 
u = r' — r and r as 

^M^ = ldp{W{u,r + u)f,ir + u,t)-W{u,r)f,{r,t)}. (62) 

Assuming that characteristic displacements are small, one can expand Eq. flB21) and arrive 
at the Fokker-Planck form of the equation for the density distribution fg{r,t), 



dfg{r,t) _ d 



A^{v)f,{v,t) + ^^{B^p{v)f,{v,t)) 



(63) 



dt dra 

The coefficients and Bap describing the acting force and diffusion, respectively, can be 
written as the functionals of the PT-function in coordinate space W. For the dimension s, 
these coefficients have the form 



A^{r) = j d'uUaW{u,r) (64) 

and 

Ba^iv) = d'uu^UpWiu, r). (65) 

For the isotropic case, the probability function depends on r and the magnitude u of the 
vector u. For a homogeneous medium, when the r-dependence of the PT-function is absent, 
the coefficient Aa = 0, while the diffusion coefficient is the constant Bap = ^apB. The 
constant B is the integral 

B = Ysj d'viU^Wiu). (66) 

This consideration cannot be applied to the specific situations in which the integral in 
Eq. fl66|) is infinite. In this case we should examine the general transport equation (I6T!) . Let 
us now consider the problem for the homogeneous and isotropic case, when the PT-function 
depends only on u. As shown in [11], after the Fourier-transformation we arrive at the 
following representation of Eq. ( 16T]) 
d/,(k,t) 



dt 



I d'u [exp(2ku) - 1] W{u)fg{k,t) ^ X(k)/,(k,t), (67) 
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where ^(k) = X{k). Let us assume the simple form of the PT-function with the power 
dependence on the distance W{p) = C/u"", where C is a constant and a > 0. Such a singular 
dependence is typical for the jump diffusion probability in heteropolymers in solutions (see, 
e.g., [26], where different applications of anomalous diffusion are considered based on the 
fractional differentiation method). For the one- dimensional case (s=l), we find 

X(k) = -4 du sin" W{u) = -23-"C| /tr-i ^C^^- (68) 

For 1 < a; < 3, this function is finite and written as 

Cr[(3-a)/2]|fcr^ 
^^^^ - "2"v^r(a/2)(a-l)' ^^^^ 



where F is the Gamma-function. At the same time, the integral in Eq. ( 16611 for such PT- 
functions is infinite, since normal diffusion is absent. 

The procedure considered for the simplest cases of the power dependence of the PT- 
function is equivalent to the equation with fractional space differentiation [10,26], 

^^^M = CA-/V,(a;,t), (70) 

Here A'^/^ is the fractional Laplacian, the linear operator whose action on the function f{x) 
in Fourier space is described by the expression A^/^/(x) = — (A;^)^/^/(A;) = —\k\^f{k). In 
the case considered above, /i = (a — 1), where < /i < 2. For more general PT-functions, 
which (for arbitrary values of the variable u) are not proportional to the power of m, the 
method described above is also applicable, although the fractional derivative method is not 
exist. 

In the case of the purely power dependence of the PT-function, the non-stationary solution 
for the density distribution describes the so-called super-diffusion (or Levy flights). The 
solution of Eq. (1701) in the Fourier space reads 

f,{k,t)=eM-C\kn), (71) 

which, in coordinate space, corresponds to the so-called symmetric Levy stable distribution 

X 



1 

(here we use the canonical notation [27]). 



_(H)Va' 



;/i,0 



(72) 
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In the general case it follows from Eq. (167|) that 

fg{k,t) = C,exp[X{k)t], (73) 

with a certain constant Ci. 

The consideration based on the PTF given in this Section allows us to avoid the fractional 
differentiation method and to consider more general physical situations of the non-power 
type of the probability transitions. Let us consider a simple example. Taking (for the 
one- dimensional case) the PT function W{u) in the form 

Win) = ci^^Ht^, (74) 
with p > 0, we arrive at the function X{k): 



oo 



Xik) = -2'~^C\kr' I d^il^^M-^i^mni^ ^ -2^--C\kr^ T{a/\k\^ a). 

Jo C" 

(75) 

It is easy to see that the function T{a/\k\P,a) is finite for 1 < a < p + 3, since for small 
distances at p > divergence is also suppressed for the powers a > 3. A simple calculation 
for a = 2 and p = 1 leads to the following result which cannot be found by the usual 
fractional differentiation method, 

T{a/\kl2) = I - arctan(|A;|/a) + ^ In [l + k'/a'] . (76) 

The asymptotic behavior of the function X{k) for k —> (or a — ?► oo) is similar, as follows 
from Eq. ( 1761) . to the case W{u) = C/u". For the case under consideration, a = 2, the limit 
X{k — )■ 0) — )• —irCk. For large values of A; (A; oo or cr 0), we find X{k) — )• {a/k) \n{k/a). 

In the general case, the universal behavior of the function X{k) is provided by the asymp- 
totical properties of the PT-function for large distances. The necessary condition for diver- 
gence of the coefficient isl<a<p + 3. 

Appendix II. Anomalous velocity diffusion for the specific case B(s)=const, 
C(s)=0 

Let us now consider formally the particular case of anomalous diffusion in velocity space, 
when the specific structure of PTF W^(q, p) provides a rapid (let say, exponential) decrease 
in the functions W'{q) and W"{q). Therefore, the exponential function under the integrals 
in the coefficients B{s) and C{s) can be expanded, which means B{s) = Bq and C{s) ^ 
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respectively. At the same time, the function W{q) = a/q" has a purely power dependence 
on q. 

In this particular case, the kinetic equation Eq. (1351) reads 



= Pos-'fg{s,t) + Bos~fg{s,t), (77) 

or, formally, in the momentum space 

= PoD'^hip,t) - Bo^[pJg{p,t)], (78) 

where u = {a — 3) (2 > z/ > 0), and we introduced the fractional differentiation operator 
D'^f{p,t) = J dss'^exp{—ips)f{s,t) in momentum space. This is a very particular case of 
the general equations Eqs. ( 122|) . ( l35i) . In this case the result is similar to the particular 
phenomenological model [17]. The stationary solution of Eq. f l77p is given by 

m = Cexp[-^] (79) 

d^sexp{—ips)exp [ — ] = / dsssin{ps)exp [ — ] (80) 

^J^o P Jo ^^0 



As an example, let us choose the case u = 1. Then, we find fg{p) 

SnC'Po 



SnC'Pp 

fM = o r/ o , ,o2/o2M2 - (81) 



In the case v = 1 the long tail of the distribution is proportional to and the distribution 
fg{p) corresponds with the Cauchy-Lorentz distribution. Normalization of the distribution 
fg{p) leads to the value C = Ug/ (27r)'^, where Ug = Ng/V is the average density of particles 
undergoing diffusion in velocity space. A similar approach can be taken for other types of 
anomalous diffusion in velocity space. 
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